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THE OH LINE CONTAMINATION OF 21 CM INTENSITY FLUCTUATION MEASUREMENTS FOR Z = 1 
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ABSTRACT 

The large-scale structure of the Universe can be mapped with unresolved intensity fluctuations of 
the 21 cm line. The power spectrum of the intensity fluctuations has been proposed as a probe of 
the baryon acoustic oscillations at low to moderate redshifts with interferometric experiments now 
under consideration. We discuss the contamination to the low-redshift 21 cm intensity power spectrum 
generated by the 18 cm OH line since the intensity fluctuations of the OH line generated at a slightly 
higher redshift contribute to the intensity fluctuations observed in an experiment. We assume the OH 
megamaser luminosity is correlated with the star formation rate, and use the simulation to estimate 
the OH signal and the spatial anisotropies. We also use a semi-analytic simulation to predict the 21 
cm power spectrum. At z = 1 to 3, we find that the OH contamination could reach 0.1 to 1% of the 
21 cm rms fluctuations at the scale of the first peak of the baryon acoustic oscillation. When z > 3 
the OH signal declines quickly, so that the contamination on the 21-cm becomes negligible at high 
redshifts. 

Subject headings: cosmology: theory— diffuse radiation 



1. INTRODUCTION 

The large-scale structure of the Universe can be ob- 
served efficiently with the intensity mapping technique, 
where the distribution of the radiation intensity of a par- 
ticular line emission from large volume cells is observed 
without attempting to resolve the individual emitters, 
galaxies, within the volume. This technique is partic- 
ularly suitable for radio observations, where the angu- 
lar resolution is relatively low. By observing the radia- 
tion at different wavelengths, the emissivity from differ- 
ent redshifts are obtained, thus revealing the three di- 
mensional matter distribution on large scales. It was 
first recognized that this method ca n be applied to the 
21cm line of the neutral hydrogen (IChang et al.l 120081 : 
iPeterson et al.1 [20091 : IChang et al.ll2010D . and provides a 
very powerful tool for precise determination of the equa- 
tion of state of dark energy by using the baryon acous- 
tic oscillation (B AO) peak of large-scale structure as 
a standard rul er ([Chang et all 120081 : lAnsari et all 120081 : 
ISeo et aO201Q[ ). More recently, it has also been proposed 
that the intensity mapping method be used for molecular 
and fine-structure lines , such as CO dGong et all [20111 : 
ICarillil 120111 ILidz et al.l I201lt IVisbal fc Loebll2010ir and 



CII (|Gong et al.ll2011bft . 

A possible problem with the intensity mapping tech- 
nique is the contamination by other lines. Unlike the 
observation of individual sources, where different lines 
from two different sources overlapping along the line of 
sight can be separated through high-resolution imaging, 
in intensity mapping the contamination from a different 
line at a different redshift cannot be easily separated. 
For optical (e.g. the Lyman a line) and for many of 
the important radio lines, there are many spectral lines 
with wavelengths longer than the line being observed, 
thus emitters at lower redshifts could become contami- 
nants and these could be obstacles in the application of 



this method. An advantage of the 21cm line for intensity 
mapping is that due to its low frequency (1420 MHz), 
there are few strong lines at a lower redshift which could 
contaminate the observations 1 . 

Nevertheless, the hydroxyl radical (OH) lines of 18 cm 
wavelength at slightly higher redshifts can potentially 
contaminate HI 21 cm observations at lower redshifts 2 . 
The A = 18 cm lines of OH correspond to four possi- 
ble transitions, with frequencies at 1612, 1665, 1667, and 
1720 MHz. Strong OH emission are produced by masers, 
originating typically in high density (n{H2) > 10 7 cm~ 3 ) 
gas near an excitation source, though the exact e nvi- 
ronm ent for the masers to happen is still not clear (|Lo I 
120051) . The 1665MHz and 1667MHz are usually much 
stronger than the other two lines, hence are named "main 
lines". In masers, the 1667MHz line is the strongest, 
whose flux is typic ally about 2 to 20 t imes greater than 
the 1665MHz line (jRandell et al.lll995l) . 

In an intensity mapping observation, the 21 cm auto- 
correlation power spectrum at a redshift z is observed. 
However, the OH emission at 1 + z' = (1 + z)^ 1 - would 
also give raise to brightness temperature fluctuations 
which can not be distinguished from the redshifted 21cm 
fluctuations in such observations. Thus, for example, the 
HI 21-cm signal at zjjj = 1,2 and 3 would be contami- 
nated by the OH 18cm emission at zoh = 1-35, 2.52, and 
3.70 respectively. Since the OH fluctuations at redshift z' 



1 See e.g. [Thompson, Moran fc Swensonl f200lT l (Table 1.1) for 
a list of important radio spectral lines, the only one below HI fre- 
quency is the 327 MHz deuterium line, which is relatively weak due 
to the low deuterium abundance. 

2 Other lines at v < 10GHz listed in the above reference 
are the CH line at 3.335 GHz, the OH line at 4.766 GHz, the 
formaldehyde(H 2 CO) line at 4.83GHz, the OH line at 6.035 GHz, 
the methanol(CH 3 OH) line at 6.668 GHz, and the 3 He line at 8.665 
GHz. These lines should be less significant than the OH 18cm line 
and we will not consider them in this work. 
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are uncorrelated with the 21 cm fluctuations at redshift 
z, the two power spectra would simply add. Although 
the OH line emission is produced with a different mech- 
anism and depends on the star formation activity, on 
large scales, we still expect the OH intensity fluctuations 
to trace the total matter densities. Its power spectrum 
should be proportional to the matter power spectrum 
at z', with a different bias factor. If not properly ac- 
counted for, this may introduce a distortion to the total 
intensity power spectrum extracted from the 21 cm ob- 
servations resulting in a shift to the BAO peaks. Given 
the low-redshift 21 cm BAO experiments are now be- 
ing developed (e.g. the proposal to conduct wider area 
surveys with GBT by building a multi-beam receiver 3 , 
the Tianlai project in China, and the CHIME project in 
Canada 4 ) it is important to estimate the magnitude of 
the potential contamination. 

For this purpose, we make use of simulations to predict 
both the 21 cm and OH intensity power spectrum from 
z = 1 to 4. We find that the contamination is generally 
small and below 1% of the rms fluctuations at z = 1 to 
3 with a large uncertainty related to the overall predic- 
tions on the OH signal. The paper is organized as follows: 
In the next section we present the method of the inten- 
sity calculation and in Section 3 we present our results 
and discuss the contamination. We will assume WMAP 
7- yea r flat ACDM cosmological model (jKomatsu et al.l 
l20ll . 

2. CALCULATION 

In order to estimate the 21 cm emission at low red- 
shi fts, we make use o f the semi-analytic simulations 
by lObrcschko w et al.l (|2009f ) . which are available as 
part of the SKA Simulated Skies 5 , and based on the 
galaxy catalog d erived from the Millennium simulation 
(|De Lucia et al.1 120071 iSpringel et "all 12005ft . These ar e 
the same simulations we have used in lGong et all (|2011f ). 
As the HI mass of each galaxy is assigned using the 
galaxy properties provided by the semi-analytical mod- 
eling of galaxy formation, we will use those neutral hy- 
drogen masses to calculate the 21 cm line intensities. 

The calculation of the 21 cm power spectr um is similar 
to what we have done in lGong et al.l f|2011h . The 21 cm 
temperature from galaxies, ass uming the signal is seen in 
emission is (jSantos et al.ll2008ft 



T b G = c(z) 
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where Ah = 0.74 is the mean hydrogen mass fraction in 
the Universe, pb = ^bPc is the baryon density and p c is 
the critical density. The c(z) takes the form as 
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The parameter c(z)xn gives the mean brightness tem- 
perature of the 21-cm emission, where xh is the mean 
neutral fraction. We assume that the neutral hydrogen 
are mostly contained within the galaxies after reioniza- 



3 https:/ /science. nrao.edu/facilities/gbt/index 

4 http://www.phas.ubc.ca/chime/ 

5 http:/ /s-cubed. physics. ox. ac.uk 



tion, the mass density pm is then given by 

Phi = / dM— M H i(M), (3) 



M min 



where dn/dM is the mass function, we take M m \ n — 
10 s M /h to be t he minimum mass for a halo to retain 
neutral hydrogen (jLoeb fc Barkanal l2001). and Af m ax = 
10 13 M /h is the maximum mass for which the gas have 
sufficient time to cool and form galaxies (the result is 
insensitive to this number). 

In above the Mm is the neutral hydrogen mass in a 
halo with mass M. The HI mass is correlated with the 
halo mass, though with some scatter. Inspired by the 
shape of the distribut ion seen in the sem i -analy tic simu- 
lation generated from [Obrcschk ow et al.l (|2009ft . we fit a 



relation of the form 



Mhi(M) = ix 1 + 
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The best-fit values of the parameters A, c±, C2, b and d 
are given in Table [TJ 

Due to the mass resolution limit of the Millennium sim- 
ulation, for M < 10 10 Mq, we cannot use the same fitting 
formula. Instead we assume Mm = x^\n b /n m )M to 

estimate the neutral hydrogen mass in halos, where Xf^ 
is the neutral hydrogen mass fraction in the galaxy. We 

set Xg 1 = 0.15 which is estimated at M = 10 10 M© 
in the simulation, and assume it does not change when 
M < 10 10 Mq. The simulation result and the best fit- 
ting curves at z = 1, z = 2 and z — 3 are shown in the 
upper panel of Fig. [1] whi ch are c onsistent with the other 
results (e.g. IMarin et al.l (|2010h ; IDuffv et al.l (pOTLft ). 

We find the HI energy density parameter fmi = Phi/ Pc 
is about 10~ 3 and insensitive to the redshift (for z < 3) 
in our calculation, which is consisten t with the observa - 
tional results (e.g IRao et al.l (|2006[ ); ILah et al.l (|2007ft : 
iNoterdaeme et all ()2009f )V Finally, we find the 21-cm 
mean brightness temperature are 481, 573 and 544 /jK 
at z=l, 2 and 3 respectively. These va lues are also consis- 
tent with an observation at z = 0.94 (Chan g et al.H2010ft 
and p revious predictions in the literature ( Chang et al.l 
12001 . 

Assuming that the 21 cm flux from galaxies is propor- 
tional to the neutral hydrogen mass Mm, then the 21 cm 
signal will follow the underlying dark matter distribution 
with a bias 



^dM^M m b{z,M) 



bm{z) 



PHI 
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where b(z, M) is the halo bias (<Sheth fc Tormenl 1 1999ft . 
The 21-cm temperature from galaxies is then T^(x) = 
T^[l + &m<5(x)], and the clustering power spectrum 
is given by P ffl (z,fc) = (f^) 2 & ffl P(fc, z). We use the 
Halof it code (S mith et al.l 12003ft to calculate the non- 
linear matter power spectrum P(k,z). Additionally, 
there is a shot noise contribution to the power spectrum 
due to the discreteness of galaxies, 
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TABLE 1 

The best-fit values of the parameters in the Mhi-M 
relation. 





2 = 1 


2 = 2 


2 = 3 


A 


2.1 x 10 s 


2.1 x 10 8 


1.7 x 10 8 


Cl 


1.0 x 10 11 


1.0 x 10 11 


1.0 x 10 11 




4.55 x 10 11 


5.60 x 10 11 


5.0 x 10 11 


6 


2.65 


2.4 


2.4 


d 


-2.64 


-2.40 


-2.25 



While the OH molecule is not as abundant as neu- 
tral hydrogen, the OH maser is brighter than the 
21 cm line intensity for the same number density of 
baryons. It is believed that OH masers are associated 
with star formation activity. OH megamasers (OHMs) 
are 10 6 times brighter than typical OH maser sources 
within the Milky Way and are found in the luminous 
infrared galaxies (LIRGs with Ltr > 10 11 Lq) and 
ultra-luminous infrared galaxies (ULIRGs with Lyr > 
10 12 L Q ) (jDarling fc Giovanellil l200l . The luminosity 
of the isotropic OHM luminosity Loh is correlated with 
the IR luminosity Llr. 

Previous studies found several relations that take the 
form as Loh oc if R , where the power in dex j3 is be- 
tween 1 and 2. For example, iBannl (|1989f) found a re- 
lation Loh <x Lf R using a sample of 18 OHM galax- 
ies. Howe ver, a flatter rela tion of £oh oc £j^ 8 was ob- 
tained by Kandalian (1996) using a sample of 49 OHM 
galaxies, after correcting for the Malmquist bias. A 
nearly linear relation between Loh and Ltr was derived 
from the Arecibo Observa tory OH megamaser survey 
(|Darling fc Giovanellil [2002] ): 



logLoH = (1.2 ± 0.1)logL m - (11.7 ± 1.2) 



(7) 



This relation has also been corrected for the Malmquist 
bias, and about one hundred OHM galaxies are used in 
this calibration. We will use this relation in our model. 

The Lir is tightly correlated with the star forma- 
tion rate (SFR) and we adopt a relat ion of the form 
(jMagnelli et al.ll20lll: iTekola et alJl25Il 



Ljr [L Q ] = 5.8 x 10 9 SFR [M yr _1 ]. 



(8) 



This relation is consistent with other works (e.g. 
lEvans et all (120061)') and have about 30% ~ 4 0% un- 
certaintv (|Kennicuttlll998l : lAretxaga et aLll2007D . 

The SFR is on a statistical sense l i nearly correlated 
with the halo mass ([Loeb et al.l 120051 : IShimasak u et al.l 
2008), with a nearly Gaussian distribution whose cen- 
tral value and variance ch anges as a function of redshift 
(jConrov fc Wec hslcr 2009). For the purpose of statistical 
calculation of OH emissivity, it is sufficient to relate the 
star forma tion rate to the halo m ass. We use the galaxy 
catalog in iDe Lucia et al.l ([2007) to derive the SFR and 
stellar mass relation SFR-M* and the star mass and 
halo mass relation M+-M as shown in the lower panel of 
Fig.[U The simulation has outputs at z = 1.38, 2.42 and 
3.57, which are fairly close to the redshifts zoh = 1.35, 
zoh = 2.52 and zoh = 3.70, which could contaminate 
the 21cm at zhi = 1, zm = 2 and zri = 3 respectively. 

We fit the SFR-M* and M+-M relations using the form 
as SFR = A x M*(l + M±/c) d and M* = B x M e respec- 
tively, and the best fit values for the parameters are listed 
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Fig. 1. — Upper panel: the Mm-M relation from the simulation 
at z = 1, z = 2 and 2 = 3. Lower panel: the SFR vs. the stellar 
mass and the stellar mass vs. the hal o mass at different red - 
shifts as derived from the galaxy catalog in Dc Lucia et al. (2007). 
The red dots are the mean value and the green dots show the ltr 
error from the simulation. The best fits are shown in blue solid 
line. 



in Table [5J For M < 10 12 M Q , analysis from the simu- 
lation indicates SFR/M ~ 10- n -10- 10 yr" 1 at 1 < z < 
4, which matches well with previous results (lLoeb et all 
120051 IShimasaku eFall l2008t IConrov fc Wechslerl I2009D 
in their applicable redshift ranges. 
We can now calculate the mean intensity of the OH 
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TABLE 2 

The best-fit values of the parameters in the SFR-M* and 
M+-M relation. 





z = 1.38 


z = 2.42 


z = 3.57 


A 


5.5 x 1CT 10 


1.2 X 1CT 9 


2.4 x 10~ 9 


c 


7 x 10 10 


9 x 10 10 


9 x 10 1(1 


d 


-1.2 


-1.4 


-1.4 


B 


5.8 x 10 2 


2.2 x 10 2 


1.2 x 10 2 


e 


0.64 


0.67 


0.69 




Fig. 2. — The OHM clustering bias and the mean Ray leigh- Jeans 
temperature at z=l, 2, 3, 4 and 5. We find the OHM signal declines 
quickly when z > 3. 



emission: 



Ion{z) = /OH 



M„ 



dM—f m (M) 47rg2 y(z)D 



A- 



'M min -*"^L 

(9) 

where the dn/dM is the halo mass function 
(jSheth fc Tormenl I1999Q . L> L and are the lumi- 
nosity distance and comoving angular diameter distance 
respectively, and y(z) = d\jdv = Aoh(1 + z) 2 /H(z), 
where x is the comoving distance, v is the observed fre- 
quency, Aoh = 18 cm is the rest-frame OH wavelength. 
The /oh is the fraction of the LIRGs together with 
ULIRGs that host the OHMs. We set fo H = 0-2 6 which 
is esti mated using the sample in IDarling fc Giovanellil 
(2002). This takes account of the fact that the OHMs 
are caused by the far-IR pumping (53 /im) in the warm 
dust (T > 45 K) which is produced an d supported by the 
star formation in LIRGs or ULIRGs (Locke tt fc Elitzurl 
2008). Note that the duty cycle does not appear in 
our formula, since it is already included in our SFR-M 
relation. 

Here we assume that the megamasers would dominate 
the contribution, so we can take M min = 10 11 M /h, 

6 This value also can be as low as 0.05, sec Klockncr (2004). 
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Fig. 3. — The rms intensity fluctuations using the power spectrum 
of the OHM and 21-cm emission. The blue dashed, dotted and 
dash-dotted lines are the OHM signal, shot-noise and total power 
spectrum. The uncertainty of the OHM emission is also shown 
in cyan region, which is estimated using the errors in Eq[7] only. 
The red lines are 21-cm power spectrum. Because the shot-noise 
power spectrum is too small compared to the signal power spectrum 
for 21-cm, the total and the signal power spectrum are almost 
overlapped together. 



as the OHMs are hosted in galaxies with molecular 
gas M H , > 4 x 10 9 (IBurdvuzha fc Vikulovl Il990l: 
iLagos et al.|[20lTl: iDuffv et al.ll201lD . The f TB (M) is the 
fraction of the LIRGs and ULIRGs for galaxies hosted by 
the halos with M > 10 11 M^/ h . whic h is estimated from 
the catalog in IDe Lucia et all (|2007D . We find f m - 1 
when M > 1.2 x 10 12 M Q , and quickly decreases for 
lower halo masses. After getting the OHM intensity, we 
can convert it into a mean Rayleigh- Jeans temperature 
Ton- 

The OHM clustering bias can be calculated in the same 
way as the HI bias (c.f. Eq. [5]), except for the weight by 
Lou instead of A/hi- The OHM bias 6oh and Ton at 
z=l, 2, 3, 4 and 5 are shown in Fig. [5J We find that the 
OHM signal declines quickly when z > 3. 

The OHM power spectrum is given by Pon(k, z) = 
I OH b QH P(k, 2), and similar to the HI case, the shot-noise 
power spectrum is given by 
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3. RESULTS AND DISCUSSION 

In Fig. |31 we plot the rms fluctuations associated with 
the power spectrum of the OH and 21-cm emission at 
different redshifts for comparison. The 21-cm signal is 
plotted in red color, while the OH in blue. In the case of 
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21 cm power spectrum, the shot noise power is relativity 
insignificant, as the number of HI galaxies is large. 

In the case of the OH power spectrum, the signal power 
spectrum evolves slowly in this redshift range, because 
the SFR is higher at high redshifts, so it counteracts the 
decrease in the matter power spectrum at high redshifts. 
The contribution of the shot noise is however very signifi- 
cant, especially at small scales (larger k), as the emission 
is mostly from the rare LIRG and ULIRG populations. 

Comparing the 21cm power to the OH power, we find 
that on the scale of the first peak of the BAO (about k = 
0.075/i/Mpc), the OH rms fluctuations are about 0.03%, 

0. 07% and 0.11% of the 21 cm rms fluctuations at zri = 

1, 2 and 3 respectively. At higher redshifts, while we do 
not show it here, we found that the 21cm signal become 
stronger as we approach the epoch of reionization, while 
the OH power become smaller and insignificant compared 
with the HI signal. 

We note that this result depends on the modeling of 
the HI and OH emission, which still has a lot of uncer- 
tainty. The exact conditions for the occurr ence of OH 
megamasers are not completely understood (|Lo I [2005). 
and the actual OH emission from a source of a certain 
star formation rate may be quite different from our model 
prediction. Moreover, the OH emissivity may not even be 
strongly correlated with the halo mass, though on very 
large scales we still expect that the OH intensity power 
spectrum to be proportional to the underlying matter 
power spectrum. The 21 cm power depends on the HI 
content of the low mass halos, which is also largely un- 
certain. 

Considering variations to our model predictions, we do 
find that in an extreme case, as shown the upper limit 
of the cyan region in Fig. 02 the OH could even super- 
sede the 21cm power spectrum, and become the major 
contribution to the observed temperature fluctuations. 
Note that we just consider the errors in the Loh-^ir 
relation (Eq. [7]) to get the uncertainty (cyan region), 
which can be greater if including the error in Ztr-SFR 



relation (Eq. HJ. Also, we note that the two relations 
above are calibrated at low redshift and tightly related 
to the redshift-dependent properties of galaxies, such as 
the galaxy metallicity. So the OH intensity may also 
increase if considering the redshift evolution effect. Of 
course, in that case, it would be more advantageous to 
use the OH emission as the tracer instead, though at 
present this does not seem to be very likely. 

As a lot of the OH power comes from shot noise, it 
may be possible to find a way to remove some of its con- 
tribution. For example, we may consider to conduct a 
targeted maser survey on the ULIRGs and LIRGs with 
sensitive telescopes which have small fields of view, and 
subtract their contributions to the temperature fluctua- 
tion. This could significantly reduce the noise power due 
to these sources. 

One way to identify and estimate the amount of pos- 
sible contamination is to cross correlate the temperature 
fluctuation at the redshift pair (24,2:2). We expect the 
temperature fluctuations should be uncorrelated, while 
contamination by OH would give raise a correlation given 
by (ST(k,zi)ST(k,z 2 )) = bm{z 2 )bou{z2)P{k, z 2 ). 

Finally, it may also be possible to make use of the mul- 
tiple lines of OH (including the four lines at 18cm and 
the lines with shorter wavelengths) to check for the con- 
tamination. If observations for individual OHM sources 
show that most of them have similar line ratios, then one 
may construct a template of OHM spectrum, and apply 
it as a matched filter on the observed spectrum to check 
for possible OH contamination. However, this would not 
be possible if observations show that the line ratios vary 
a lot. 

We thank Prof. Xingwu Zheng for helpful discus- 
sion. This work was supported by NSF CAREER AST- 
0645427 at UCI, by the 973 program No. 2007CB815401, 
the NSFC grant No. 11073024, and the John Temple- 
ton Foundation at NAOC; and MGS and MBS ac- 
knowledges support from FCT Portugal under grant 
PTDC/FIS/100170/2008. 



REFERENCES 



Ansari, R. et al., 2008, arXiv:0807.3614 
Aretxaga, I. et al., 2007, MNRAS, 379, 1571 
Baan, W. A. 1989, ApJ, 338, 804 

Burdyuzha, V.V. & Vikulov, K.A. 1990, MNRAS, 244, 86 
Carilli, C, 2011, ApJ, 730, 30. 

Chang, T., Pen, U., Peterson, J. B., McDonald, P., 2008, 

Phys. Rev. Lett., 100, 091303. 
Chang, T., Pen, U., Bandura, K., Peterson, J. B., 2010, Nature, 

466, 463. 

Conroy, C. & Wechsler, R. H. 2009, ApJ, 696, 620 
Darling, J. & Giovanelli, R. 2002, AJ, 124, 100 
De Lucia, G. et al., 2007, MNRAS, 375, 2 

Duffy, A. R., Kay, S. T. , Battye, R. A., B ooth, C. M.,Vecchia, C. 

D. & Schaye, J. 2011. larXiv: 1107.37201 
Evans, A. S., Solomon, P. M., Tacconi, L. J. & Downes, D. 2006, 

ApJ, 132, 2398 
Gong, Y., et al. 2011, ApJ, 728, L46 
Gong, Y., et al. 2011b, arX iv: 1107.35531 
Kandalian, R. A. 1996, 39, 237 
Kennicutt, R. C. 1998, ARAA, 36,189 

Klockner, H.-R. 2004, PhD Thesis, University of Groningen 
Komatsu, E., et al. 2011, ApJS, 192,18 

Lagos, C. P., Baugh, C. M., Lacey, C. G., Benson, A. J. Kim, H. 

& Power, C. 2011, larXiv:1105. 2294 
Lah, P., et al. 2007, MNRAS, 376, 1357 
Lidz, A., et al. 2011, arXiv: 1104:4800 
Lo, K. Y., 2005, ARA&A, 43, 625. 



Lockett, P. & Elitzur, M. 2008, ApJ, 677, 985 

Loeb, A., & Barkana, R. 2001, ARA&A, 39, 19 

Loeb, A., Barkana, R. &c Hernquist, L. 200 5, ApJ, 620, 553 

Magnelli, B., et al. 2011. larXiv:1101. 24671 

Marin, F. A., Gnedin, N. Y., Seo, H. & Vallinotto, A. 2010, ApJ, 
718, 972 

Noterdaeme, P., Petitjean, P., Lcdoux, C. & Srianand, R. 2009, 

A&A, 505, 1087 
Obreschkow D., et al. 2009, ApJ, 703, 1890 
Peterson, J. B. et al., 2009, arXiv:0902.3091 
Randell, J., Field, D., Jones, K. N., Yates, J. A. & Gray, M. D. 

1995, A&A, 300, 659 
Rao, S. M., Turnshek, D. A. & Nestor, D. B. 2006, ApJ, 636, 610 
Santos, M. G. et al., 2008, ApJ, 689, 1 
Seo, H. J. et al., 2010, ApJ, 721, 164. 
Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 
Shimasaku, K., Hayashi, M., Motohara, K., et al. 2008, ASP 

Conf. Ser., 399, 75 
Smith, R. E., et al. 2003, MNRAS, 341, 1311 
Springel, V., et al. 2005, Nature, 435, 629 

Tekola, A. G., Vaisanen, P., & Berlind, A. 2011. larXiv:1101.3495l 
Thompson, A. R., Moran, J. M., Swenson, G. W. Jr., 2001, 

Interferometry and Synthesis in Radio Astronomy, John Wiley 

& Sons Inc., New York. 
Visbal, E. & Loeb, A. 2010, JCAP, 11, 016 



